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Abstract. It is often the case in mathematical analysis that solving an open 
problem can be facilitated by finding a new set of coordinates which may il- 
lumniate the known difficulties. In this article, we illustrate how to derive an 
assortment coordinates in which to represent the relativistic Boltzmann colli- 
sion operator. We show the equivalence between some known representations 
I15II31| . and others which seem to be new. One of these representations has 
been used recently to solve several open problems in [~~ 
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Dedicated to the memory of Naoufel ben Abdallah and Carlo Cercignani 

1. Introduction and main result 
The relativistic Boltzmann equation can be written as 

p»d„F = C(F,F). 
In this expression the collision operator [3j[T5] is given by 

C(f,h)= ( ^1 %l %W(p,q\p',q')[f(p')h(q')-f(p)h(q)}. 
Jr™ Q Jr™ Q je» P 

The transition rate, W(p, q\p' , q'), from [3] in n-dimensions (n > 2) is denoted 

I / ft \ 3— n 

(1) W(p,q\p',q') = -^) S *(Q,e)6 il+n Hp» + q »-p^-q / n, 

where <j{q, 9) is the differential cross-section which is a measure of the interactions 
between particles. This is an important model for fast moving particles. Standard 
references in relativistic Kinetic theory include [T3lll5, 28.40,41]. In this paper we 
give a complete reduction of the collision integrals for the operator C(f, h), deriving 
several sets of coordinates for the particles momentum, some old and some new. 
The rest of our notation is given after the following historical discussion. 



2000 Mathematics Subject Classification. Primary: 76P05; Secondary: 83A05. 

Key words and phrases. Relativity, Boltzmann, collisional Kinetic Theory, Kinetic Theory. 

R.M.S. was partially supported by the NSF grant DMS-0901463. 



1 



2 



R. M. STRAIN 



1.1. Historical discussion. The relativistic Boltzmann equation is the primary 
model in relativistic collisional Kinetic theory. In the next few paragraphs, we will 
provide a short review of the mathematical theory of this equation. We mention a 
few books on relativistic Kinetic theory as for instance [13l[T5l[2"8U40[|4Tj . 

Carlo Ccrcignani and collaborators wrote several research works on relativity 
and the Boltzmann equation including [6ll7ll9l-fT4]: many of these works were in 
collaboration with Gilberto Medeiros Kremer [9UT3]. Together [9] they showed 
that the original proof of the formula for the "summational invariants" , due to 
Boltzmann, also applies in the special relativistic case. In addition Cercignani- 
Kremer wrote a book on the relativistic Boltzmann equation |13j . This follows the 
book by Cercignani-Illner-Pulvirenti |8j on the Newtonian Boltzmann equation. 
The author learned about the physics of the relativistic Boltzmann equation from 
the books of Cercignani-Kremer [13] and de Groot-van Leeuwen-van Weert [15] . 

Now early results include those of Lichnerowicz and Marrot |37j . who wrote 
down the full relativistic Boltzmann equation, including collisional effects, in 1940. 
Later, in 1967, Bichteler [4] showed that the relativistic Boltzmann equation has 
a local solution. Then Dudyhski and Ekiel-Jezewska [20], in 1988, proved that 
the linearized equation admits unique solutions in L 2 . Afterwards, Dudyhski |19] 
studied the long time and small- mean- free-path limits of these solutions. 

In the context of large data global in time weak solutions, the theory of DiPerna- 
Lions |17| renormalized solutions was extended to the relativistic Boltzmann equa- 
tion in 1992 by Dudyhski and Ekiel-Jezewska [21] . This result uses the causality of 
the relativistic Boltzmann equation [22][23]. Results on the regularity of the gain 
term are given in [U|47]; the strong L 1 compactness is studied by Andreasson pQ. 
These are generalizations of Lions [38j result in the non-relativistic case. Further 
developments on renormalized weak solutions can be found in [35ll36j . 

Let us point out studies of the Newtonian limit [5 ]l44j for the Boltzmann equation. 
We mention theories of unique global in time solutions with initial data that is "near 
Vacuum" as in [24 l l29 ] l44] . Andreasson, Calogcro and Illner [2] prove that there can 
be blow-up in the presence of only the gain term. Notice the study of the collision 
map and the pre-post collisional change of variables from [30]. Then [25] provides 
uniform instability estimates for the relativistic Boltzmann equation. Now there 
is a mathematically rigorous result connecting the relativistic Euler equations to 
the relativistic Boltzmann equation via the Hilbert expansion as in [39] . 

Next we discuss several results for initial data that is nearby the relativistic 
Maxwellian. In 1993, Glasscy and Strauss [21] first proved the global existence and 
uniqueness of smooth solutions on the torus T 3 . They also established exponential 
convergence to Maxwellian. Their assumptions on the differential cross-section, 
cr(g, 8), covered the case of hard potentials. In 1995, they extended that result to 
the whole space case [32j where the convergence rate is polynomial. Subsequent 
results with reduced restrictions on the cross-section were proven in [25] , using the 
energy method from for instance [33]; these results also apply to the hard potentials. 
Recently, the author [45] has developed a weighted L°°(T 3 x R 3 ) theory of global in 
time unique solutions for the soft potentials, and rapid convergence to Maxwellian, 
under the general physical assumption proposed in [181120] . 

We point out results on global existence of unique smooth solutions which are 
initially close to the relativistic Maxwellian for the relativistic Landau-Maxwell 
system [42], and then for the relativistic Landau [27] equation as well. Further 
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[49] proves the smoothing effects for relativistic Landau-Maxwell system. And [48] 
proves time decay rates in the whole space for the relativistic Boltzmann equation 
(with hard potentials) and the relativistic Landau equation as well. 

1.2. Notation. In this section we define several notations which will be used 
throughout the article. A relativistic particle has momentum p = (p 1 , . . . ,p n ) £ R™, 
with its energy defined by p° = y/c 2 + \p\ 2 where \p\ 2 = p ■ p. Here c denotes the 

speed of light. We use the standard Euclidean dot product: p ■ q = f Y^i=\P l 1 l - 
As is customary we write p M = (p°,p) where p M also denotes the //-th element of 
(p°,p). In general, Latin (spatial) indices etc., take values in {1, . . . , n}, while 
Greek indices k, A,/i, is, etc., take on the values {0,1,..., n}. Indices are raised 
and lowered as usual with the Minkowski metric g M „ and its inverse g^ v , where 

(<?£ii>) == diag(— 1 1 • • • 1) is an (1 + ?i) x (1 + n) matrix. In other words p^ = g^, v p v '. 
We furthermore use the Einstein convention of implicit summation over repeated 
indices. The Lorentz inner product is then given by 

Notice p^p^ = — c 2 . Now 9 M = (c~ 1 dt, V x ) so the streaming term of the relativistic 
Boltzmann equation is p^d^ = c podt + p ■ V x , where V x is the spatial gradient. 
Conservation of momentum and energy for elastic collisions is expressed as 

(2) p^ + qf 1 =p'^ + q^. 

These conservation laws are enforced by the 1 + n delta functions in ([T]) . The angle 
9 in the Boltzmann collision operator (JIJ is then defined by 

(3) cos0 =(p^- ^(p^-^Ve 3 . 

Note that this angle is well defined under ©, see [551 Lemma 3.15.3]. Here the 
relative momentum, g = Q^p^yq 11 ), is denoted 

(4) q d = f yV-^XjV-^) = ^2(-p^ At _ c 2) = ^(pOgO-p-g-c 2 ). 
Furthermore the quantity s = sijp^^q^ 1 ) is defined as 

(5) S d = f -(pf + q»){p, + q^) = 2 {-p»q, t + c 2 ) > 0. 

This is in other words s = 2(p°q° —p- q + c 2 ). Notice that s = g 2 +4c 2 . To proceed, 
we will quickly review the Lorentz transformations. 

A digression on Lorentz transformations. In this subsection, we discuss a few 
elementary aspects of Lorentz transformations which will be useful throughout the 
rest of this paper. Let A be a (1 + n) x (1 + n) matrix (of real numbers) denoted by 

A = (A M !/ )o<p, 1 /<n- 

For the basics of Lorentz transformations, we refer to [13] and [46] . 
Definition 1. A is a (proper) Lorentz transformation if det(A) = 1 and 

g K \ A A „ = g^, (/i,i> = 0,l,...,ra). 
This implies the following invariance: p K q K = p K gn\q X = (A K M p /i )<7 K A(A A , y 9 I/ ). 
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We will use the notation A to exclusively denote a Lorentz transformation. Any 
such A is invertible and the inverse matrix is denoted A -1 = (A^)o<n,u<n so that 
A" re A K = 5" ' , where S v is the standard Kronecker delta which is unity when the 
indices are equal and zero otherwise. It follows from Definition [T] that 

(6) A; = A\ g Kll . 

Definition [T] further implies that (A _1 ) M 1V = A J 1 is a Lorentz transformation. 

In this paper, we are exclusively concerned with (proper) Lorentz transforma- 
tions A, depending on p M and which further satisfy 

(7) AW + q v ) = {V~s^,...,Q). 

We recall that s is defined in Note that there are several Lorentz transforma- 
tions which will satisfy ([7]). In particular, in dimension n = 3, three such Lorentz 
transformations were given in [441 Appendix A] . 

The condition (J7J means that you are mapping the particle momentum to the 
"center of momentum" frame p + q = (this is sometimes called alternatively the 
center of mass frame). Notice that (J7J further implies 

(8) A° M (p^-^)=0. 

This follows quickly from the fact that (0 implies A 1 jp^ 1 = —A 1 ^ for any i <E 
{1, . . . ,n}. Then Definition [Hallows us to further deduce (JSJ as a result of the mass 
shell condition: -c 2 = p K p K = (A K ll p fJ ')g K x(A x J/ p v ). 
We give a standard example of these mappings. 

The Boost matrix. The most common Lorentz transformation is probably the Boost 
matrix. Given v = (v 1 , . . . , v n ) £ R n , the (1 + n) x (1 + n) Boost matrix is 

def / p -pv 

-pv l n + ( p _ 1 W 



where p = (1 — |v| 2 ) ^ 2 and l n is the n x n identity matrix. Notice that A^ has 
only n free parameters. Our goal is to choose v such that A& satisfies ([7]). Let 

dot P + q def p° + q° 

p v + <r v s 

Plugging these choices into At above, we obtain that 

/ P°+q" p+q \ 

^ Afc = I S (p+q)®(y+q) )■ 

\ ~ i« + - ij " i; +g v y 

By a direct calculation this example satisfies (|7|)- 

1.3. Main results. We will now state our main results in the language of the 
notation just introduced. For this we use the operator Q{f, h) = f C(f, h)^ as 

(10) Q(f,h) = ^f §/ %[ %W(p,q\p',q')[f(p')h(q')-f(p)h(q)}. 

The main point is to carry out the reduction of the number of integrations in this 
expression by evaluating the 1 + n delta functions from ([T]). In the literature, 
there are two approaches to performing this goal [HJ[2T]. One [3T] uses algebraic 
manipulation of polynomials; this results in a representation (Theorem U) that is 
well known in the mathematics literature and has been widely used (at least in 
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dimension n = 3). The other representation |15j (from Corollary [3] and Theorem 
[2} , which starts by using the change of variables as in , is sketched in physics 
texts but seems to be hard to locate in the mathematics literature on the relativis- 
tic Boltzmann equation; this approach was developed (as in Corollary [3] below in 
dimension n = 3) in the authors thesis [43]. Additionally the representations of 
the collision operator (fTU| from Corollary [3] have only recently been used to solve 
several open problems on the relativistic Boltzmann equation 34,39,44,45 . 

For this reason, we are motivated at this time to write down a complete math- 
ematical proof of these different representations. Our results improve upon those 
given previously [15,31 in the following ways. We prove these representations in n 
dimensions with n > 2. Note that the representation in Theorem|4]is dimensionally 
dependent; there is an additional term fBJfl which is not present when n = 3. Fur- 
thermore it is discussed in the physics paper [3] that there are physical situations in 
which the Boltzmann equation may be of interest in dimensions other than n = 3. 
We also prove the exact formula for the post-collisional energies as in (|2"Tj) . We give 
the precise expressions for the angles in (fl~3|) . ([M]) . and (|22j) which do not seem to 
have been previously computed. What the author finds most interesting is that 
we can show that there are several alternative expressions for the post-collisional 
momentum and energy, as in (|15p and (|16p . one for each Lorentz transformation 
satisfying (|7|) . We hope that the availability of these additional alternative represen- 
tations may be useful to future investigations in the relativistic Boltzmann theory; 
in particular we observe in Corollary [5] the equivalence of these different represen- 
tations. This equivalence has been crucial to our recent proof (joint with Yan Guo) 
of the global in time stability for the relativistic Vlasov-Maxwell-Boltzmann system 
[34| with near Maxwcllian initial conditions. 

To begin we state the center of momentum |7J reduction. 

Theorem 2. (Center of momentum reduction). Recall (|10p and |T]). For any 

suitable integrable function G : l n X R" X 1" X 1" -» R, it holds that 

(n) Jzf %\ %w(p, q \p', q ')G( P , q , P >, q >) 



dw v a(g,0) G(p, q ,p, q ), 
where lu = (to 1 , . . . ,U) n ) <G § n_1 and v = v {p,q) is the M0ller velocity given by 

I 12 ) v = v {p, q ) = 7-5-0. 

4 p u q" 



The angle 9 in this expression is defined by (|13|) and (jl4[) . The post-collisional 
momentum and energy above are defined by (|15|) and (|16j) respectively. 



The angle in the reduced expression in Theorem [2] is defined by 

k 

(13) cos9 = — • ui, 

\k\ 

where k G R™ is given as 

(14) k id ^A i ^- q n=A i (p -q )+A i j (pi- q i), (i = l,...,n), 
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Moreover the post-collisional momentum satisfy (i = 1, . . . , n) 
(15) 2 2 

Here we use the tensor L^ K = f — g^ x A K X . Furthermore, the energies are 



(16) 



x 

= L oo^_ L (« Wj | > 



10 rOO , T 0j, , 6 

q =L — + L J uj j -. 



2 J 2 

As usual, in each of these expressions, we implicitly sum over j G {1, . . . , n}. Note 
that the formulae above hold for any Lorentz transformation A satisfying only ([TJ . 
Now H]), flTU|) and (JTTJ) together imply that 

Q(f,h) = [ dq [ duv a(g,9) [f(p')h(q') - f(p)h(q)}, 



where the relevant quantities are defined above as in (flU|) . (fT4")l , (fT5|) and (|TB)) . 
If we use the Lorentz boost © satisfying ([7]) we have the following simplification. 

Corollary 3. In the particular case of the Boost matrix the post-collisional 
momentum (|15p in the expressions above can be written precisely as 

, p + q . q ( . , , u , . (p + g) • u 

/ p + q q ( , / 1W , ^(p + ?)' w 
9 = — -H w + (p_ 1)(p + 9) TT^ 

where p = (p° + q°)/^/s. Furthermore, the post-collisional energies are given by 

p o + q o Q 
P =^- + ^{ P + q)-^ 

, p° + q a g 
q = ^i--2^ {P + q) ' UJ - 
Additionally in the angle (|13|) . the vector (|14p can 6e simplified to 

, P+JZ, o On , / \ I / 1V , s (g + g) • (g - g) 

fc = -- "^-(P -?) + (p-9) + (p- 1)(P + 9) ^p^l • 

In i/iese formula we use s from (|5|) . 

Now we turn to the expression given by Glassey-Strauss in |31j . 
Theorem 4. (Glassey-Strauss reduction). Alternatively 

P 'w{p,q\p\q')G{p,q,p\q') 

duj ^cfo^ B n{p,q,u) G(p,q,p',q'), 

where B n (p,q,uj) is given by (|19[) and (p',q') on the right are given by (|20l) . TTie 
ang/e # is also defined by 



p 9 jr" q jr™ p 
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In the above reduction, we consider the expression 

f (pP + gO) W | w .^_^ 



(18) B(p,q,u;) d =c 



Then the kernel in Theorem |4] is given by 

(19) B n (p, q, u) = B(p, q, u) (A(p, q, W ))"~ 3 . 

where 

ef (P° + q>°q° |w • - 



[(p0+g0)a_( w .[p +g ])2] 

Notice that B^(p, q,u)) = B(p,q,uj) when the dimension is n = 3. 

For the reduction in Theorem [4j the post-collisional momentum are given by 

p' =p + a(p,q,uj)uj, 
q = q- a(p,q,uj)uj, 

where 

2(/+g>.(A-g"p)} 
a[p,q,u)- {p o + q0) 2 _ {uj . [p + q]) 2 ■ 

And the energies can be expressed as p'° = p° + N° and q'° = q° — N° with 

dc f 2u ■ (p + q){u) ■ (p°q - q°p)} 



(21) N l 



(pO + g o )2 _ (w . b + g])2 



These formula clearly satisfy the collisional conservations ©. The angle © in 
Theorem 21 can then be reduced to 

(22) cos, = l- 8 



e 2 (p + g °) 2 -( W -b + g]) 2 



Evidently, in this case is not simply a dot product of a unit vector with oj. 
Moreover, assuming the collisions are elastic as in ([3]), we have the invariance: 

u-{p q-q p)=U>-(q> p>-p'°q'). 

Therefore, for fixed uj € S" _1 , B n (p, q, cj) = B n (p', q 1 , u). Then this kernel is invari- 
ant under pre-post change of variables (p 1 ,q') — > (p,q) as in [30] . 
In this case (p} , (|10j) and Theorem 0] together illustrate that 

Q(f,h) = [ dq [ duj S -^§2p- B n (p,q,u>) [f(p')h(q') - f(p)h(q)], 

where the quantities inside the integral are defined as in (f2"2"j) . (f^Uj) . (fH?)) and (|T5|) . 
We mention now the following useful corollary. 



Corollary 5. Combining Theorem^ and Theorem^ yield, 



■s 



duj v <t(q,6) G(p,q,p',q') = [ duj S<7 ^' ^ B n (p,q,Lo) G(p,q,p' ,q'), 

where the angle and post-collisional momentum on the left are defined as in Theorem 
[JJ and the similar expressions on the right are given as in Theorem^ 
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Now in the variables (|20|) . it was computed in [30l Theorem 1] that the mapping 
(p,q) — > {p',q') has the following Jacobian 

d(p',q') 



d(p, q) 



Polo 

poqo' 



Note also that when (p, q) —> (p' , q') then additionally (p' , q') — > (p, q). We conclude 

dp [ dq [ du S -^-§ 2 f- B n (p,q,u) G(p,q,p',q') 

dp dq duj — r^r— B n (p,q,(j) G(p' ,q' ,p,q). 

J Rn J Sn -i p°q° 

This holds with the variables (|20[) which are used in Theorem [4j However iterating 
this formula and using Corollary [5] we additionally observe that 

(23) dp dq du v <r(g,9) G(p,q,p',q) 

Jk" Js." Js"- 1 



dp J dq J doj v a(g,6) G(p',q',p,q). 

The variables in these integrals are those from (|15j) which are used in Theorem [2] 
This article is organized as follows. In the next Section [5] we will prove Theorem 
[21 from which we conclude Corollary [3J Then in Section [3] we will prove Theorem 
SI We remark that Corollary [5] will be an important part of |34j . 

2. Center of Momentum reduction of the Collision Integrals 

In this section we will prove Theorem [2] and in particular . To this end we 
consider the following integral 

( 24 ) 7 =/ %[ %s<?(Q,0)5( 1+n \p» + q»-p'»-q'nG( P ^q^p'^q'n 

ir q jr™ p 

= [ %[ %MQ,0)S (1+n) (P l * + <l fl -p'> 1 -<fnG(p,q,p',q'), 

jr™ q jr" p 

where 9 is defined by (J3)). Above and in the following, for convenience, we write 
G(p^,q fJ, ,p'^,q' 11 ) = G(p,q,p' ,q') when there is no opportunity for confusion. 
Our goal will be to prove that 

(25) I = 2 2 -" / dcu g n - 2 ^S a{g, 9) G(p, q,p', q'). 

JS"- 1 

Above we use the definitions fig]) , (fill) . fI5"]) . and ([IB]). Notice then, by ((!} and 
(|T2"j) . that this implies ([TTI) . 



Integral reduction. We now reduce the collision integrals in the center of mo- 
mentum system. We first claim that / from (|24p can be written as 



I= l %l %MQ,o)6 {1+n \^(p" + ql-p^-q'nG(p^q»,P , ^q , n- 
jr" q jr™ p 

Here A'^ should satisfy ([7]). In the expression above from © we have 
(26) p'» d =* Afp'" = A\g KV p' v , and = g>* A\ g KV q w . 
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This holds because ^ and are Lorentz invariant measures and 

(AM (p" + q v - p ,v -q' v )) = 8 {1+n) {p» + - p^ - g"*). 
But notice that the claim is not true unless the angle 9 from ([3]) is redefined as 

cosO = {A^(p» q»)W v q' v )/Q 2 = || • 

where fc is defined in (| 14|) . We have also employed the following calculation 

= yV - g")0v - %) = y/WUp« - q K )}g^Wxip x - q x )} = \k\. 

For that we used ((HJ). Note also, as a result of ([7]), we can deduce (from the delta 
function in /) that p' + q' = which further yields p'° = q'°. 

Then the integration over q' can be carried out immediately and we obtain 

1 f dp' 



I 



2 Un (p'O) 2 

Here now, with p' = —q', the angle is 

k p' 
cos 8 = 2—— ■ —. 

\k\ g 

Next change to polar coordinates as p' = \p'\us with us S S™ -1 and 

dp' = \p'\ n - l d\p'\djjj. 
We use the following calculation to compute the delta function 



2 2 V iyi+v/c 2 + |p'| 5 

|g 2 -H 2 

We have just employed s = + 4c 2 from Q and ([5]). Thus 

i r _{\q-\p'\){\q + \A) 

V s P 



2 iVi+v^+Fp 



Now if = then s = g 2 + Ac 2 grants that y/c 2 + \p'\ 2 = \\fs. Therefore, 



We plug this in to see that 

1 /-oo /. I nn-1 3/2 /i \ 

/= 2 i, ^Q,e)S^-Q-\p'\jG{p^q^p^,q'^ 

i \ n f n—1 3/2 

- / dto 7 l——a(g,e)G(p^ q ^pi^ q in 
2/ Js—i (iVs) 2 £> 



d W Q n - i s x l i <x(Q,6)G<j?,q!>,it' x t cli x ). 
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Now consider w M = (p'°, \p'\u>) = ^(y/s, Qw) and tD" = (p'°, — \p'\oS) = \{\fs, —Qui). 
Then in the integral for / above, we replace ([26]) with 

(27) p'» = g» x A\ g KV Q" , q'» = g» x A\ g KU Q» . 

We calculate below that these are (|15|) and ([T6]) . Furthermore the angle is given as 
in $33]). We thus obtain JUJ, (HI]) and Theorem [2 

Post-collisional momentum. We now calculate the post-collisional momentum 
and their associated energies (|15|) and (p~6|) . We use (f27| to compute 

p'" = -L^g^uT = -L^Q K = L»°^- - L^|, 

where i MK = f -g' lA A K A . Similarly 

h 2 J 2 

These establish that ([27]) is equal to ([UJ and (JTSJ). 

Notice that, using (l27l) . we recover the identities of conservation for elastic col- 
lisions ([2]) from (j6j) and ((Jj as follows 

+ g") = 0, ... ,0) = cD" + £}" = A^p'" + 

The last equality follows from ([2~T]) since 

AW + q' u ) = K g uX A K A g Ka {uj a + u a ) = uj» + w", 

which itself is implied by A^ v g uX A K X = g^ K ; the previous statement is a consequence 
of Definition [1] To obtain ^ we apply the inverse, as in ([5]), to observe 

p>» + J* = A»A K v (p' v + q'n = A»A K v (tf + tt) =P» + <?. 

This completes the proof of Theorem [2] 

3. Glassey-Strauss reduction of the Collision Integrals 

The goal of this section is to prove Theorem @] This reduction was given by 
Glassey & Strauss |31| in dimension n — 3 and without the presence of the ar- 
guments, (g,6), in the differential cross section, a(g,9). The reduction below is 
essentially similar to |31j in that we perform the examination of the roots of poly- 
nomials; it is also however slightly different from |31| in that we integrate the radial 
variable, r, below over the entire real line R (rather than r > 0). We also observe 
a dimensionally dependent factor A(p, q, ui) as in (|19p which is not present when 
n = 3. We are further able to compute the formula for the angle 9 in (|22]) and 
establish the equivalence of the representations as in Corollary [5] 

We will use the expression for I in ([24]) . We claim that 

(28) / = - c (|)"" 3 i dco sa(g, 9) B{p, q, lo) (A(p, q, co)) n ~ 3 G(p, q,p',q'), 

where (p',q') are defined above as in ([20]) . In this case, using ([24]) . (Q]), ([T8]) and 
(fig]) , we obtain Theorem H] subject to (J2SJ). 

To reduce the number of integrals in (f2~4"]) . we split I = |/ + ^1. Letting q' = 
p + q — p' we can immediately remove n of the delta functions in the first copy of I. 
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Next translate p' — > p + p' so that g' — > q — p'. Then switch to spherical coordinates 
as p' = rw, <ip' = r n ~ 1 drdu where r £ [0, oo) and u 6 S™ _1 . We obtain 

/ —£^sa( g ,e)S(p' + q ' -p - q °)G(p, q ,p', q '), 

where p' = p + ru and q' = q — ru. 

Now consider the other \l. Let p' = p + q — q ' to remove n of the delta functions. 
Next translate q' — > g + q' so that p' — > p — g'. And switch to spherical coordinates 
as g' = rw, <ig' = r n ~ l drdu where r <G [0, oo) and u E § n_1 . Then 

k I= k / —^sa( g ,e)S(p' + q' -p - q )G(p, q ,p',q'), 
L A Jo JS"- 1 P 1 

where p' = p — ru and q' = q + ru. Further change variables r — > — r so that 
\ I= \f i M "Jot^ sa{ Q ,9) 6(p> + q> -p°- q )G(p,q >P >, q >), 

where now p' = p + ru and g' = g — rw. 

We combine the last two splittings to conclude that 

i = \i+\i= \ I** I W ~So^ Me,0) % /0 +g /0 -P -g°)g(P,g,P / ,g / ) 1 

where p' = p + ru and g' = q — ru. Now the angle ([3]) satisfies 
_ {p o_ q o )(p/0 _ q/ o ) + ip _ q) . {p ,_ q , ) 



(29) cos 6» = 



<? 2 



(p° - g°)(p'° - g'°) + \p- q\ 2 + 2r(p -q)-u 
Q 2 



We will return to this expression below. 

We now focus on the argument of the delta function. For Aj > (i = 1, 2), we 
use the identity S(Xi — A2) = 2Ai<5(A 2 — A|) and ([2]) to sec that 

S(p'° + q'° -p°- q°) = 2(p° + q Q )6((p'° + g' ) 2 - (p° + g ) 2 ) 
= 2(p° + g°) <5(2p'V - {(p° + g ) 2 - (p'°) 2 - (g'°) 2 }). 

If {(p° + g ) 2 - (p'°) 2 - (g'°) 2 } > then this is 

= 8(p° + g°)p'°g'° <5(4(p'°) 2 (g' ) 2 - {(p° + g ) 2 - (p'°) 2 - (g'°) 2 } 2 ). 

If {(p° + g ) 2 — (p'°) 2 — (g'°) 2 } < 0, then the delta function is zero. Now we write 
the argument of the last delta function above as 

P(r) = -(p° + g ) 4 + 2(p° + g°) 2 {(p' ) 2 + (g'°) 2 } 

-{(p'°) 2 + (g' ) 2 } 2 +4(p'°) 2 (g' ) 2 

= -(P° + q a ) 4 + 2(p° + g°) 2 {(p' ) 2 + (g'°) 2 } - {(p'°) 2 - (g'°) 2 } 2 . 
Plugging in p' = p + ru and g' = g — ru we observe that 

(p'°) 2 - (g'°) 2 = |p + rw| 2 - |g - ru\ 2 = (p ) 2 - (g ) 2 + 2ru • (p + g). 
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This means that p(r) is quadratic in r. Moreover, 

P(0) = -(P° + I") 4 + 2(P° + q°) 2 {( P Q ) 2 + (q Q ) 2 } {(P ) 2 (q Q ) 2 } 2 

= -(P° + I*) 4 + 2(P° + q°f{(p°) 2 + (q°) 2 } (P° + ?°) V ^ q°f 
= -(p° + <Z°) 4 + (P° + <?°) W) 2 + (<Z°) 2 } + 2p°g (p° + g ) 2 = 0. 

We conclude that p(r) = AD\r 2 — 8D 2 r for some Di,D 2 £ R. 
We will now determine D\, D 2 . Expanding 

(P'°) 2 + {q'°? = 2c 2 + \p + ruj\ 2 + \q- ru\ 2 = (p ) 2 + (q ) 2 + 2ruj ■ (p - q) + 2r 2 . 

We plug the last few calculations into p(r) to write it out in terms of r as 

p{r) = -(p° + q ) 4 + 2(p° + q°) 2 {(p ) 2 + {q ) 2 + 2ru ■ (p - q) + 2r 2 } 
-{(p ) 2 -(q°) 2 + 2ru-(p + q)} 2 
= -(p° + q°) 4 + 2(/ + g°) 2 {(p°) 2 + (q°) 2 + 2ru • (p - q) + 2r 2 } 

- {(P°) 2 ~ (<7°) 2 } 2 - ^ 2 {u • (p + q)} 2 - 4r{u ■ (p + q)}{(p ) 2 - (q ) 2 }. 
Rearrainging the terms 

p(r)=4{(p + qy-{u-(p + q)} 2 }r 2 

+ 4 {(p° + q°f{u ■ (p - q)} -{oj-(p + q)}{(p ) 2 - (q ) 2 }} r 
- (p° + qy + 2(p° + q°) 2 {(p°) 2 + (q°) 2 } - {(p ) 2 - ( 9 °) 2 } 2 
= 4{(p° + qy-{o J -(p + q)} 2 }r 2 

+ 4 {(p° + qy{u ■ (p - q)} -{u J -(p + q)}{(p ) 2 - (q ) 2 }} r. 
Equivalcntly we have 

£i = {(/ + <7°) 2 - Wp + <?)} 2 }, 
2D 2 = -(p° + q ) 2 ^ ■ (p - q)} + {u-{p + q)}{{py - (q ) 2 }. 

We further calculate D 2 as 

2D 2 = -{{ P y + (qy + 2p°q°){u ■ (p - q)} + {u-(p + q)}{( P y - (q ) 2 } 
= w ■ p{-2(q ) 2 - 2p°q°} + lj ■ q{2( P y + 2p°q°} 
= 2{(p° + q °)u J -(p q-q p)} 

Thus, p(r) = 4Dir 2 — 8D 2 r with these definitions. 

We plug this formulation for p(r) into the full integral to obtain 

1 f +oa f 

I=- / \r\ n - l drdu s(t{q,0) 8(p° + g°)<5(4D 1 r 2 - 8D 2 r)G(p,q,p' ,q'). 

Then we extract the factor 4£>i from the delta function to obtain 

1= / {rr'drdu s<j(q,9) [P ^ Q >- 6(r{r - 2D 2 /D 1 })G(p,q,p' ,q'). 

J-oo JS"- 1 U l 
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We use the identity 8{r{r - 2D 2 /Dx}) = 
delta function drops out because of the |r 



2D 2 
n-1 



{S(r) + S(r - 2D 2 /Dx)}. The first 
factor (n > 2). We thus have 



/ = / °° / Irl^drdw ^(g, 0) ^ +^ } <5(r - 2D 2 /Dx)G(p, g,p' , g')- 
We have used Di > 2 [301 : note Z?2 can be negative. Evaluate the delta function 



Notice that from (jTT?)) and (TT5)) we have 

2"" 2 ( (P° + 9°) = § (f ) " 3 9, «) (^(p, 9j -))' 1 - 3 . 
This grants us (|28p. and thereby Theorem 2) Furthermore, the angle (|2"5)) becomes 

-(p° - 9°)(p'° - + \p- q\ 2 + Hp - 1) • w 



cos 6* 



-(P° - g°) {(p '°fcy )2} + ^ - ^ + 4f (P - 9) 



_(pP _ gO) {(/) 2 -(y+2 r .(p + ,)} + [p _ g |2 + 4g?(p _ g) . 

P 2 

Above we have used ((2]) and r = 2D 2 /Dx- Further calculations yield 

-(P° - 9°) 2 - fej 4§J W • (p + g) + |p - g| 2 + 4^(p - g) ■ w 

(30) 008 0= 1^—^ 

g 2 

D 2 u-(p + q) + (p-q)-LJ 

2 



Di Q 

=i+4^r 2_ 



Dx g 2 

1 AuJ-{p\~q*p)Y 



Dx 6 2 



Note dnOJ) is exactly 



Post-collisional energy. Wc show in this part how to compute the post-collisional 
energy ([/H]) . The point is to write p'° = p° + N° and q'° = q° — N°, and also recall 
from p'° — c 2 + \p'\ 2 and q'° = y/c 2 + \q'\ 2 that one can obtain 

-c 2 = p'% = -(p° + N ) 2 + \p'\ 2 = -(p° + N ) 2 + \p + a{p,q,u J )u\ 2 . 

Similarly, 

-<? = q'% = -(g° - iV°) 2 + |g - a(p, q,u>)u\ 2 . 
By expanding both of these expressions, and subtracting the result, one obtains 

2(p° + g°)iV - 2u ■ (p + q)a(p, q, u) = 0. 

Solving this equality for N° establishes (j2"Tj) . A related computation (but for a 
different purpose) can be found in Cercignani-Kremer [T3l Section 1.4.3]. Q.E.D. 
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Remark. Wc point out that there is unfortunately a misprint in our recent paper 
[44] if the dimension n > 2 is not n = 3. In particular the transition rate W at the 
top of the paper [33], should be replaced by the transition rate from ([T]) (and [3]). 
The main difference between the two is the factor g 3 ~ n , which is unity when n = 3. 
Furthermore, the expression (1.9) in [44] for the collision operator is correct when 
n = N = 3 (N is the notation for the dimension used in [44]) , but otherwise the 
kernel B(p,q,ui) in [44] (1.9)] needs to be replaced by B n (p,q,uj) in ([19]) . In other 
words the factor (A(p, q, u))) n ~ from (fT9]l is missing from [44] (1.9)]. We point 
out that this factor in the expression [44j (1.9)] does not affect the main theorems 
of [33] ■ Furthermore the condition in [331 (2.7)] on the collisional cross section is 
written as < 7 < -3 (which is empty); this [331 (2.7)] should be < 7 < N. 
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